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In this paper, we study the inequality indices for some models of wealth exchange. We calculated 
Gini index and newly introduced fc-index and compare the results with reported empirical data 
available for different countries. We have found lower and upper bounds for the indices and discuss 
the efficiencies of the models. Some exact analytical calculations are given for a few cases. We also 
exactly compute the quantities for Gamma and double Gamma distributions. 
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I. INTRODUCTION 


Socio-economic inequality GHi is manifested in the existence of unequal rewards and opportunities for social posi¬ 
tions or statuses in a society. Structured, recurrent patterns of unequal distributions of goods, wealth, opportunities, 
rewards and punishments are mainly measured in terms of inequality of conditions, and inequality of opportunities. The 
former refers to the unequal distribution of income, wealth and material goods, while the latter refers to the unequal 
distribution of ‘life chances’ of individuals. This is somehow reflected in measures such as level of education, health 
status, and treatment by the criminal justice system. Socio-economic inequality often results in crisis, political unrest 
and instability, conflict, war, criminal activity and finally affects economic growth Q. Initially, economic inequalities 
were studied in the context of income and wealth p-y, but the notions and observations have led to widespread 
research, see e.g. Ref. for various socio-economic inequalities. The study of inequality in society is a 

topic of global focus and utmost current interest, bringing together researchers from various disciplines. 

By the end of the 19th century, Pareto 0 made extensive studies and found that wealth distribution in Europe 
follows a power law for the rich, commonly known to be the Pareto law. Subsequent studies have revealed that the 
distributions of income and wealth possess some globally robust features (see, e.g., 0): the bulk of both the income 
and wealth distributions seem to reasonably fit both the log-normal and the Gamma distributions. Economists have 
a preference for the log-normal distribution [l^, [l^ , while statisticians 0 and physicists root for the 

Gamma distribution for the probability density or Gibbs/exponential distribution for the corresponding cumulative 
distribution. The high end of the distribution, known as the ‘tail’, is well described by a power law as observed by 
Pareto. Formally, the probability distribution of wealth is given by 


P[m) 


F{m) 

m.^+" 


for m < me, 
for m > TOc, 


( 1 ) 


where a is a constant and v is called the Pareto exponent, ranging between 1 and 3 0 (See Ref. for a historical 
account of Pareto’s data and some recent sources). F{m) is some function which could be exponential. Gamma or 
lognormal. The crossover point me is extracted from the numerical fittings. 

One of the key class of models uses the kinetic theory of gases [2l| , where the gas molecules colliding and exchanging 
energy was mapped to agents meeting to exchange wealth, following certain rules [l9j |. In these models, a pair of 
agents agree to trade, each save a fraction A of their instantaneous money/wealth and exchanges a random fraction 
of the rest at each trading step. The distribution of wealth in the steady state, Pirn) matches well with the empirical 
data. When the saving fraction A is fixed, i.e., in case of homogeneous agents (CC model hereafter) [2^, P{m) are 
very well approximated to Gamma distributions [^ . It is important to note that, in reality, the richest follow a 
different dynamic where heterogeneity plays the key role. To obtain the power law distribution of wealth for the 
richest, one needs simply to consider each agent as different in terms of the fraction of wealth he/she saves in each 
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FIG. 1: (a) Gini and fc-index schematic; the solid red line is the Lorez curve, the cumulative proportion X of ordered individuals 
(from lowest to highest) in terms of the cumulative proportion of their sizes Y. Gini index is given by the ratio g = A/{A-\- B). 
The fc-index is computed from the intersection (fc, 1 — fc) of the reverse diagonal with the Lorenz curve, (b) Gini index from 
World Bank data for several countries over years. 


trading [^ . which is very natural to assume, because it is quite likely that agents in a market think differently from 
one another. With this very little modification, one can explain the whole range of wealth distribution [I^. When A 
is distributed uniformly in [0,1) and quenched, (CCM model hereafter), i.e., for heterogeneous agents, one obtains a 
Pareto law for the probability density of wealth P{m) m~'^ with exponent v = 2 [i3)[2J|. Several variants of these 
models, find possible applications in a variety of trading processes (n 1^. 

Socio-economic inequalities are quantified in various ways. The most popular measures are absolute, in terms of 
indices, e.g., Gini [I^, Theil [2^, Pietra and the recently introduced k index [2^. The alternative approach is 
a relative measure, in terms of probability distributions of various quantities, but the most of the above mentioned 
indices can be computed from the distributions. Most quantities often display broad distributions, usually lognormals, 
power-laws or their combinations. For example, the distribution of income is usually an exponential followed by a 
power law [2^ (see Ref.Q for other examples). 

To compute the Gini index, one has to consider the Lorenz curve [s^, that represents the cumulative proportion 
X of ordered individuals (from lowest to highest) in terms of the cumulative proportion of their sizes Y (See Fig. [TJi). 
X can represent income or wealth of individuals. The Gini index {g) is defined as the ratio between the area enclosed 
between the Lorenz curve and the equality line, to that below the equality line. If the area between (i) the equality 
line and the Lorenz curve is A, and (ii) that below the Lorenz curve as B, the Gini index is given hy g = A/(A + B). 
The recently introduced ‘k index’ is defined as the fraction k such that {1 — k) fraction of individuals possess k 
fraction of income or wealth (See Fig.[TJi) [sij. 

In this paper, we investigate the inequality in wealth in some models of wealth distribution which are inspired by 
kinetic theory of gases. We mainly discuss the results for two well studied models (CC and CCM) and a new model 
for bimodal distribution of wealth. We numerically compute the inequality indices, Gini index and A:-index to quantify 
the inequalities. Gini index g, the most popular and widely used measure for inequality in case of income and wealth 
distribution, can take value from 0 to 1. The value g = 0 refers to complete equality and 5 = 1 represents completely 
inequality. The meaning of /c-index, say, for a wealth distribution, is the following: k fraction of the top wealthiest 
people possess 1 — k fraction of total wealth. We found that in both CC and CCM models there are some upper 
and lower limits of the indices. For CC model, g varies between 0 and 0.5, and k from 0.5 to 0.68. Similarly, for 
CCM model, g varies between 0.4 to 0.85. Therefore, both models independently do not cover the possible theoretical 
range of values of q and k. We find that the range of g as found from empirical data (0.2 — 0.7) (see Fig. [TJ using 
World Bank data [^) can be well covered by CCM model. We also considered a model where two groups of agents 
have fixed but different saving propensities. Depending on the combinations, the resulting probability distribution 
of wealth is found to be unimodal or bimodal. The phase boundaries, depending on the ratio of the two groups and 
the combination of values of their saving propensities are also computed numerically. The bimodal distribution seems 
to fit well to a combination of two gamma distributions (double -Gamma distribution). Gini index and k index are 
calculated for this model for different combination of parameters. Next, we considered gamma and double-Gamma 
distribution and computed certain quantities like Lorenz curve and Gini indices. 
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II. MODELS AND NUMERICAL SIMULATION RESULTS 

Kinetic exchange models of wealth distributions [l^ serve as simple paradigmatic models for exchange of wealth 
in an economy. The main idea is that agents possess wealth which is redistributed upon trading with others. The 
‘economy’ is assumed to be a ‘closed’ one, in the sense that neither the number of agents change nor does that total 
amount of wealth in the system, and the economic activity is limited to exchange of wealth according to certain rules. 
The basic model in the framework is just the random sharing of wealth, motivated by random exchange of energy 
between gas molecules, as in the framework of kinetic theory of gases Q- The basic money exchange model 
imitates the kinetic exchange in an ideal gas, but subsequently developed models incorporate the notion of ‘savings’. 
In the following, we discuss these models, and also compute the inequality measures like Gini and fc-index. 


A. CC model 

Savings come as a natural ingredient in a trading economy. In each trading step, a pair of agents exchange their 
wealth in the following way: they keep a fixed fraction A of their wealth to themselves and the rest 1 — A fraction 
is pooled up to be randomly split among the two [1^. In this model (CC model hereafter), agents are homogeneous 
- all of them save the same fraction of their instantaneous wealth at each trading step. Formally, the dynamics is 
defined by 


mi{t + 1) = Xrriiit) + r(I - A) + mj{t)) 

mj{t + 1) = + (I - r)(l - X){mi{t) + 

where r is a random fraction in [0,1], drawn in each time (exchange) step. mi{t) and mi(t + I) are the wealth 
of the fth agent at trading times t and (t + I) respectively. The ‘saving propensity’ A is a fixed fraction in [0,1). 
A = 0 corresponds to complete random exchange (DY model) while A = 1 gives no dynamics. For A = 0, Pirn) = 
exp(—m/(m)) is exponential, for which Lorenz curve, Gini and k index were derived [Hj. However, for finite A, 
P{m) = Cm°‘ exp(—m/T), has a form of Gamma distribution [2^, where T = and C = —• The exponent 

a is related to the parameter A as a = We plot P(m) vs. m for different values of A in Fig [2^. For these 

simulations (and for each case discussed in the paper), the average wealth (m) is set to unity. We measured inequality 
in the distributions in terms of Gini and k index and plotted in Fig. |5]3 for different values of A. 




FIG. 2: (a) Probability distributions P{m) for different values of A in the CC model, (b) Gini and fc-index for the entire range 
of A in CC model. Simulation results are shown for N = 1000 agents with average wealth unity. 
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FIG. 3: Probability distributions P{m) in the CCM model for various distributions of A, given by n(A) = (1 + (5)(1 — A)'^. The 
power law exponents are v = 2 + S. (a) For <5 = 0,0.5, 1.0, 2.0; (b) For negative values of <5. (c) Gini and fc-index for a range of 
J values in the CGM model. Simulation results are shown for N = 1000 agents with average wealth unity. 


B. CCM model 

In this model [H, [111 (CCM model hereafter), agents are assumed to be heterogeneous, the saving fraction A for 
each agent is different, drawn from a given distribution 

n(A) = (1 + ^)(1 — A)^, with — 1 < (5 < oo, (3) 

where A is a fraction in the interval [0,1 — e], where e is arbitrarily small and positive. The dynamics of exchange 
follows 

mi{t + 1) = + r((l - Xi)mi[t) + (1 - Xj)mj{t)) 

mj{t+ 1) = Xjmj{t) + (1 — r)((l — Xi){mi{t) + (1 — Xj)mj{t)), 

where r is a random fraction in [0,1], drawn in each time (exchange) step. Xi is the saving fraction of agent i whose 
wealth is mi(t) at trading step t. Xi are quenched and drawn randomly from n(A) (Eq. [^. The asymptotic form of 
the steady state distribution of wealth is given by [l^ 

P{m) ^ (5) 

We compute P(m) vs. m for different values of 5. In the simulations, we allow values of A until a certain fixed upper 
cutoff (0.999 for our case; i.e., e = 0.001), because the population of agents with A close to 1 is already high, and we 
have to restrict any agent assuming a saving propensity very close to unity. As a result, there is no finite size effect 
in the calculated values. 
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In Fig. |3^ we plot the distributions P(m) for non-negative values of 5. We find that the probability distribution 
Pirn) follows Eq. [S]for most of the range of m, until an exponential cut-off, which is a result of the truncation of n(A) 
close to 1. In Fig.|3lD we plot the same for negative values of 5. The probability distribution P{m) again follows Eq. [5] 
for most of the range of m. However, it is quite interesting to note that for 5 < —0.5, P(m) shows a second peak at a 
large value of m, say m*. This m* moves towards (m) = 1 as <5 —>■ — 1. This is quite easy to explain theoretically: as 
(5 —> — 1, n(A) is peaked near A —^ 1, essentially more and more fraction of agents have very high saving propensities, 
a situation similar to A I in CC model. We recall that this situation will tend to produce P{m) peaked at average 
money per agent ((m) = I), which is equivalent to more “equality”. P{m) = A(m — (to)) with (to) = 1 for A = I in 
CC model; here A(-) is the Dirac delta function. Then Gini index g = 0 and k = 0.5. In comparison, A are distributed 
in CCM model and making <5 —— 1 makes n(A) further peaked near A = 1, more agents have ‘similar’ values of 
saving, close to unity and produce the second peak close to to = (m) = 1. The peak moves towards (to) as S ——1. 
Additionally, Gini index g ^ 0. It may be noted that, because of this modification over the power law (Eq. [S|) in the 
distribution function P(m), the standard relationship between Gini index and the Pareto exponent (see e.g.. Ref. 0) 
is not valid here. 

We measured inequality in the distributions in terms of Gini and fc-indices and plotted in Fig. [31(c) for different 
values of S. Inequality seems to be maximum around S « —0.5 


C. Model for bimodal distribution and phase diagram 




FIG. 4: (a) Wealth distribution P{m) for two groups of agents of same size, Ni = N 2 = 500, i.e., p = 0.5 with saving 
propensities Ai = 0.2 and A 2 = 0.8. The total system size is Ni + N 2 = N = 10®. The steady state wealth distribution is 
fitted to a combination of two Gamma distributions + a 2 m°‘‘^e~^^'^ with parameter values ai = 0.98, 02 = 10.57, 

Pi = 4.82 and P 2 = 7.23. (b) Phase boundaries for various values of p, showing regions of unimodal and bimodal distributions. 
The simulations have been performed for N = 10® agents. 


Bimodal distributions in wealth distributions are not uncommon and is also observed in firm sizes [s^. In 
the following we propose a very simple modification in the kinetic exchange model framework, to produce bimodal 
distribution of wealth. 

Let us now consider two groups of Ni and N 2 agents, with saving propensities Ai and A 2 respectively. Let p = 
iVi/(iVi -h A 2 ). Agents’ saving propensities remain unchanged over time, and they exchange money using the same 
rule as CC model (Eq. |2|). For example, let us consider the case p = 0.5 i.e., iVi = A 2 . Let Ai = 0.2 and A 2 = 0.8. 
After exchanging their money, the steady state distribution is shown in Fig.|4l The distribution clearly shows bimodal 
distribution. We fit the distribution with a combination of two Gamma distributions aiTO“'^e“^i™ -I- a2TO“^e“^’^™ 
with parameter values ai = 0.98, 02 = 10.57, /3i = 4.82 and P 2 = 7.23 as shown in the Fig 0^. We found 
that all combinations of (Ai,A 2 ) do not give the bimodal distributions. The combination of values giving bimodal 
distribution are shown in the phase diagram. The boundary region is roughly estimated for various values of p 
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(= 0.1,0.2,0.3,0.4,0.5) and shown in Fig. Idjo. We also computed Gini index and fc-index for different combination of 
(Ai, A 2 ) for various values of p (Fig. [5]). 



FIG. 5: Top: Gini indices for different cases where N\ agents have saving fraction Ai and N 2 agents have A 2 each, for different 
values of p = Ni/{Ni + A^ 2 ). Bottom: fc-index for the same. The simulations have been performed for Ni + N 2 — N = 10^ 
agents. 


III. GAMMA DISTRIBUTION AND ITS INEQUALITY STATISTICS 


For the CC model, the steady state wealth distribution closely fits gamma distributions [23|. Let us compute the 
inequality measures considering such a distribution. For the Gamma distribution: 


P(m) oc m“e /3 = 1/T 

we evaluate the inequality statistics. The cumulative distribution is given as 

x{r) = [ m°^e-^^dm = 

Jo 

where "f{a,x) is an incomplete gamma function defined by 

7 ( 0 , a:) = f 

Jo 

Hence, we have the normalization constant of the distribution ([61) as 

r(a + l) 


c(oo) = 


pl+0 


where we define the Gamma function by 


pOC 

r(a) = / t°'~^e~*dt = j{a,oo). 

Jo 


( 6 ) 


(7) 


( 8 ) 


(9) 


(10) 
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FIG. 6 : Inequality in Gamma distributions: (a) Lorentz curve for several values of a; (b) Gini index 3 as a function of a. 


Thus, we have 

X( I = = 7 ( 0 : + l,/3r) 

x{go) r(Q: + 1 ) 


( 11 ) 


Similarly, we have 


and 


Hence, 


y{r) = 


f/3r 


/ 32+0 




7 (a + 2,/3r) 


2/(00) 


r(a + 2 ) 

^2+a 


Y( 1 = 7(« + 2,/3r) 

- y{oo) r(a + 2) 


( 12 ) 


(13) 


(14) 


Therefore, the Lorentz curve is given by m and dm) . In Fig. [B^, we plot the Lorentz curve for several values of 
a. We also checked numerically that the curve is independent of the parameter /3. Next, to calculate Gini index, we 
check if the following relation is satisfied: 


7(0 + I, a:) = 07 ( 0 , x) — x°‘e~^ 

Applying it to our case, we immediately obtain 

7 (a + 2 , jSr) = (a + 1 ) 7(0 + 1 , fir) — {fir)°‘^^ e~ . 


This reads 


r(a + 2 )y = (a + l)r(a + 1)A - {I3r)‘^+^ q-. 


If we notice r(Q; + 2) = (a + l)r(a + 1), we obtain 


X-Y = 


(/3r)«+ie-/5’- 

r(a + 2 ) 


Accompanying the derivative 


dX _ (I3r)‘^e-^^ 

dr r(a: + 1) 


(15) 

(16) 
(17) 


(18) 


( 19 ) 

















with X — Y, we get the Gini index g as 


= 2 


= 2/3 


[ {X-Y)dX 

Jo 

L 


(^^)a+le-/ 3 r 


0 r(a + 2 ) 

2/3 


r(a + i) 


-dr 


r(a + 2)r(a + l) Jo 


poo 

Jo 


r(2(a + l)) 


22“+ir(a + l)r(a + 2)’ 


( 20 ) 


which is independent of (3 and we recover the exponential case by setting a = 0 as 

r( 2 ) 1 

2 r(i)r( 2 ) 2 ' 

Note r(l) = 1. The cases of a = 1, 2 are given by 

r(4) 3 

“ 23r(2)r(3) “ 8 
^ 7 ( 6 ) ^ 

“ 25r(3)r(4) “ 16 


( 21 ) 


( 22 ) 

(23) 


where we used r(a + 1) = ar(a) recursively. In Fig. [SI we plot the g as a function of a, and values of g corresponding 
to a = 1, 2 are also indicated. 


IV. MIXTURE OF GAMMA DISTRIBUTIONS: UNIMODAL AND BIMODAL DISTRIBUTION 



FIG. 7: (a) Mixture of Gamma distributions as a bimodal distribution lEg. I24II . (b) Lorentz curve for the same (from Eg. 1261 
and [27]). 


We next consider the mixture of Gamma distribution as 


P{x) = (I - /)x"^ 


Kl —1 




r(Ki) 0 r 




/*2 —1 






(24) 


Obviously, for / = 0, the unimodal Gamma distribution takes its maximum at a; = 0 i(ki — I). However, for a 1, 
the peaks are located at the x as solutions of 


1 -/ 


_ (0i(«i - I) - x) (02(a^2 - 1) - x) x«-ie-"/^^ = 0 


nKi)ei 


n^2)0^2 


(25) 
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Hence, the locations of the peaks are dependent on the choice of parameters 0i, 02, ^ 2 ,/■ For the Lorentz curve, 

we obtain 


where we defined 


Xir) = + f^9'(^2,?’/02) 

+ rI^^27(K2 + 1, r/ 02 ) 
~ Y{oo) 


Y (oo) 


(l-/)0i 


r(/ci + 1 ) 
r(«i) 


+ /02 


r(«:2 + i)l 
r(« 2 ) J- 


(26) 

(27) 


(28) 


In Fig [Tl^, we plot the Lorentz curve for several choices of the parameters. For the above Lorentz curve, the Gini 
index g(/, ki,/ t 2 ,0i, 02 ) is calculated as 


■[ d7’?''='"^e"’'/®l{7(K2,r/02) - T7^7('«2 + 1,r/02)} 

Jo F(oo) 

pOO Q 

/ drr'"=“^e“’'/®={7(Ki,r/0i) - / . 7(^1 + l,r/0i)} 

Jo F(oo) 

^ -1 pOO Q 

) ^ Jo - y^7('«2 + l,r/ 02 )}. 


0 ^^r(Ki)r(«; 2 ) 
/(!-/) 
e--r{^,)T{K2) Jo 

f V 1 

r(K 2 ) 


(29) 


Knowing the parameters /, ki, K 2 , 0i, 02, one can compute g numerically from the above expression. 


V. DISCUSSION 

Empirical data [s^ shows that Gini index varies mostly in 0.2 — 0.7 (see Fig.[TlD). The GC model gives Gini index 
in the range 0 — 0.5. The inequality decreases monotonically with increasing saving propensity X. g = 0.5 for A = 0, 
the wealth distribution P{m) is a perfect exponential distribution, giving the maximum value of inequality for this 
model. On the other extreme, when A —> 1, P{m) approaches a Dirac J-function A(m — (m)), for which g ^ 0 
(Fig. [2|). Hence this model does not reproduce most of the range of real Gini indices. In reality, Gini index rarely go 
below 0.3, but often goes beyond 0.5. 

In GGM model, however, the range of the Gini index is quite wide, and in fact, overlaps with almost the entire 
range of empirically observed Gini index values. In fact, in the asymptotic limit of Eq. [5]for practically infinite value 
of (5, P(m) should approach an uniform distribution in [0,1], which will yield a value of Gini index equal to 1/3. 
In Eig. [3 }d, we observe that for large values of <5, there is a tendency to saturate to a value close to 0.4, which we 
anticipate, might as well approach 1/3 for S —^ 00 . 

In some empirical analysis, the bulk of the wealth distribution resembles Gamma distribution. We analytically 
computed the Lorenz curve and the Gini index for Gamma distributions. There are even some instances where the 
wealth distribution are found to be double peaked [ 2 ^. We propose a variation of the kinetic exchange models to 
model this, and a combination of Gamma distributions to fit the resulting distribution. The steady state distribution 
is unimodal or bimodal depending on the combination of values of the saving propensities and the relative fraction of 
agents of the two groups. The phase boundaries for specific cases have been computed using numerical simulations. 
We also show that it is possible to derive an exact expression for the Gini index, considering Gamma distribution as 
the best fit to the numerically computed steady state wealth distributions. 

The critical studies of kinetic exchange models of wealth distributions seem to yield more and more interesting 
aspects, not only in terms of theoretical understanding of the models, but also when com par ed to empirical data. 
For instance, one of the recent studies explain city size statistics using the same framework |34l |. There has not been 
only a very few studies (ssj that discuss inequality measures in reference to models. Further research will be able to 
elucidate the usefulness of such a simple framework in understanding complex socio-economic phenomena. 
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